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Abstract 

We consider a classical analogue of the well known quantum two- 
slit experiment. Charged particles are scattered on flat screen with 
two slits and hit the second screen. We show that the probability 
distribution on the second screen when both slits are open is not given 
by the sum of distributions for each slit separately, but has an extra 
interference term that is given with the quantum rule of the addition of 
probabilistic alternatives. We show that the proposed classical model 
has a context dependence and could be adequately described with 
contextual formalism. 



1 Introduction 



It is well known that the classical rule for the addition of probabilistic alter- 
natives 

V = V!+V 2 (1) 

does not work in experiments with elementary particles. Instead of this rule, 
we have to use quantum rule 

V = V 1 +V 2 + 2<fP{P 2 cos6. (2) 

The classical rule for the addition of probabilistic alternatives is perturbed 
by so called interference term. The difference between 'classical' and 'quan- 
tum' rules was (and is) the source of permanent discussions as well as various 
misunderstandings, see e.g. on general references [1]- [19] . We just note that 
the appearance of the interference term was the source of the wave- viewpoint 
to the theory of elementary particles; at least the notion of superposition of 
quantum states was proposed as an attempt to explain the appearance of a 
new probabilistic calculus in the two slit experiment, see, for example, Dirac's 
book [1] on historical analysis of the origin of quantum formalism. We also 
mention that Feynman interpreted (|2]) as the evidence of the violation of the 
additivity postulate for 'quantum probabilities', [5]. 

In particular, this induced the viewpoint that there are some special 
'quantum' probabilities that differ essentially from ordinary 'classical' prob- 
abilities. We also remark that the orthodox Copenhagen interpretation of 
quantum formalism is just an attempt to explain (Q) without to apply to 
mysterious 'quantum probabilities'. To escape the use of a new probabilis- 
tic calculus, we could suppose that, e.g. electron participating in the two 
slit experiment is in the superposition of passing through both slits. We 
mention that, in particular, this implies that quantum particles do not have 
trajectories. 

However, there is another approach to quantum experiments that is not 
so strongly based on special "non-classical" features of elementary particles. 
This is so called contextualist approach. In experiments with elementary 
particles we have to take into account whole experimental arrangement, see 
N. Bohr [3] and W. Heisenberg [4]. Thus quantum probabilities are context- 
depending probabilities. Here the term context is used for a complex of 
experimental physical conditions. The contextualist approach to quantum 



2 



mechanics was developed in many directions, see e.g. [9] -[19]. Recently the 
classical probabilistic derivation of quantum rule (0) was presented in the se- 
ries of papers of one of the author's [20] -[23]. This derivation demonstrated 
that it seems that special quantum features are not important to get inter- 
ference modification (^) of classical rule ([!]). Interference can be induced for 
macroscopic systems by variations of context. It is not important what kind 
of physical systems, micro or macro, are prepared by a complex of physical 
conditions. Theoretical investigations [20]- [23] demonstrate that we could, 
in principle, get interference for macroscopic systems. 

These theoretical considerations stimulated the numerical investigation 
presented in this paper. We consider a classical analogue of a well known 
quantum two-slit experiment. Charged particles are scattered on flat screen 
with two slits and hit the second screen. We show that the probability dis- 
tribution on the second screen when both slits are open is not given by the 
sum of distributions for each slit separately, but has an extra interference 
term that is given with the quantum rule of the addition of probabilistic 
alternatives. In principle, we can introduce complex amplitudes of (classi- 
cal!) probabilities and work with macroscopic quantities in the Hilbert space 
framework, cf. [22]. 

2 The Model 

We consider a classical analog of the two slit experiment (Fig.l). The uni- 
formly charged round particles are emitted at point e with fixed velocity with 
the angles evenly distributed in the range [0, 2tt). Each particle interacts with 
the uniformly charged flat screen Si . The charge distribution on the particle 
and the screen stays unchanged even if the particle comes close to the screen. 
Physically this is a good approximation when the particle and the screen are 
both made of dielectric. There are two rectangular slits in the screen (on the 
Fig.l the slits are perpendicular to the plane of the picture). Particles pass 
through the slits in screen Si and gather on screen S2. 

We consider three experiments. In the first one the bottom slit is closed 
with the shutter, in the second - the upper slit, and in the third both slits 
are left open. The charge distribution on the shutter is the same as on the 
screen, i.e. in the first two experiments one can think as if the uniformly 
charged screen has only one slit. In this and several paragraphs below by the 
screen we mean screen Si. 



3 



2R 




21 


x 

e 






D 


d 



Figure 1: Two-slit experiment. 

Now let us write the equations of motion in each of three experiments 
(2 = 1,2,3) 

mr = Fi (3) 

where r determines place of the particle. Here F^ is force affecting the 
particle in each experiment. It is given by the Coulomb's law 

^=Swv\ ds (4) 

where r' is a vector from an element on the screen to the particle, q is charge 
of the particle, a is charge density on the screen, i.e. charge of a unit square. 
We integrate over the surface of the screen, the integration region T>i is plane 
of the screen except the splits, as it was mentioned above it is different in 
each experiment depending on which slits are opened. 

Projecting equations ©-(§]) to xy-plane, where x and y denotes horizontal 



4 



and vertical coordinates of the particle respectively we get 
mi = qa J dy' J dz' 



x 2 + (y- y') 2 + z' 2 ) 3 / 2 

f f y-y' ^ 

my = qa J dy' j dz - 



{x 2 + (y- y') 2 + z' 2 fl 2 

where Tj indicates the integration region for the i-th experiment. In our 
previous notations V>i = Tj x IR. We have 

r 1 = (-00,/) u (/ + 2^+00) 

r 2 = (-00, -l - 2R) U (-1, +00) (6) 
T 3 = (-00, -/ - 2/2) U (-/, l) U (Z + 2i2, +00) 

Here 2Z is the distance between slits and 2R is the height of the slit. 
Integrating the rhs of (H) we get 



mx = qa 

{a,b)cTi 



E( b — y a — y\ 
2 arctan arctan 
\ x x J 



my = qa In (x 2 + {b — y) 2 ) — In (x 2 + (a — y) 2 ) 

(a,b)cTi 



(7) 



where the notation (a, b) C Tj means that the sum extends over all subranges 
of Tj given in (^j. For example for i = 1 we have two summands with 
(a = — 00, b = I) and (a = Z + 2i?, 6 = +00), and (0) will take the following 
form 

.. n ( l-V l + 2R-y 

mx = 2qa it + arctan arctan 



x x 



my = qa In 



x 2 + (l-y) 2 (8) 



x 2 + (I + 2R - y) 2 



here we took into account that arctan(±oo) = ±7r/2 and the sum of the 
logarithms for a = —00 and b = +00 vanishes. 
We take the following initial values 

x(0) = — D x(0)=f cosa 
?/(0) = ?/(0) = sin a 



5 




where angle a is a random variable uniformly distributed in [0, 2tt). The con- 
stant parameters Vq and D are initial velocity and distance between emitter 
and the screen. 

Particles are emitted at point e (see Fig.l), move obeying (0),(0) passing 
through slit(s) in the screen Si and gather on the screen 5*2. Having points 
where particles hit the screen S2 we compute frequencies with which particles 
appear on screen S2 as a function of coordinates on the screen. We interpret 
this frequencies as probability distributions. We are interested in computing 
the probability distribution over a vertical line on screen S2 with z = 0. That 
is why we consider a motion only in the xy-p\a.ne and initial values (|9[) do 
not contain ^-coordinate. 

We solve the equations of motion (0) with initial conditions (f|) numeri- 
cally. We use Rungie-Kutta 4th order switching to Adams 4th order method. 
We used GNU C++ (g++) compiler to realize the simulation on Ultra- 
SPARC computer running Solaris. We had to explore about 10 5 trajectories 
and we used 4-processor parallel computer located an Vaxjo University. The 
computation process was easy to make parallel as moving particles are not in- 
teracting, i.e. one could think as if they were emitted with long intervals. The 
algorithm automatically adjusted the computation precision making shorter 
steps when the particle comes near to the first screen or the coordinates 
(x and y) changed more than minimum precision allowed. The first stage 
of computation was calibration when the algorithm determined the angle 
ranges for which the particles passed through the slits and hit the second 
screen. This reduced the angle range from [0, 2tt) to a set of ranges, which 
are different in each experiment. In fact we used symmetry of the first two 
ones (when only upper or lower slit is opened) making computations only for 
the first one. The second screen was separated with cells of equal size, the 
diameter of a particle. The number of particles which hit into each cell was 
calculated and interpreted as a probability distribution. 
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Figure 3: Interference term, cos6 = (2P i2 — (Pi + Pq)) / y/P~iP~2- 

Let us denote the probability distribution in the first experiment (only 
upper slit is opened) as Pi = Pi(y), in the the second experiment (only 
lower slit is opened) as P 2 = P2O/), and in the third experiment (both slits 
are opened) as P i2 = Pi 2 (y). Although since the force is different in each 
experiment, see (]3|), it is quite clear that (Fig. 2) 

*Wt + T (10) 

To become an equality the above equation should have an extra term 

Pi2 = y + y + V / ^cos0 (11) 



where \fP\P2 cos 6 is a so-called interference term (Fig.3), and 6 = 6(y) is 
spread along y-axis. 
The function 

2Pi 2 -(Pi + P 2 ) 
cosfl = . 12) 

is shown on the (Fig.2). Please note that as there are ranges where Pi or 
P 2 are equal to zero, i.e. P1P2 = the function cos 9 is not determined and 
from ( |TTD we see that P i2 does not depend on it. 

Conclusion. We have shown that the proposed classical model has a 
context dependence and could be adequately described with contextual for- 
malism. Quantum like behavior for macro systems is demonstrated. We sim- 
ulate quantum-like interference for macroscopic objects. Such a simulation 
essentially reduced the gap between micro and macro worlds. [] 

In particular, recent experiments of the group of A. Zollinger [23] and the Boulder- 
group [24] can be interpreted as successful steps in this direction. 
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